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Abstract 

We consider a one dimensional model of an electron in a doubly (or nearly) degenerate band that 
interacts with elastic distortions. We show that the electron equations of motion reduce to a set 
of coupled non-linear Schrodinger equations. For the case of interband electron-phonon coupling 
stemming from local Jahn- Teller interactions, new multicomponent self-localized polaron solutions 
-vector polarons- are described and classified. The phase diagram for the different types of vector 
polarons in this model is presented. By interpreting the components of the orbital doublet as those 
of spin-i, our results can also be used to describe bound magnetic polarons. 
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Nonlinear excitations are known to play an important role in shaping the optical and 
transport properties in quasi-one-dimensional electron systems. The best-studied ex- 
ample is the kink soliton in trans-polyacetylene, as elucidated in the seminal work of Su, 
Schrieffer, and Heeger P]. Many studies have concluded that there are a number of tech- 
nologically intriguing materials whose properties are determined by solitons, polarons and 
jipolarons such as conducting polymers 11,0,01, organic Q and high-Tc superconductors 
J, halogen-bridged mixed valence transition metal (MX) chains ja,0i, and the mangan- 
ites y, Recently discovered nanoscale materials such as polymeric fuUerides and carbon 
nanotubes may well be added to this list. 

We consider the effects of local Jahn- Teller (E^bi) and pseudo-Jahn- Teller ((a4i©i?i)®6i) 
electron-phonon coupling in a continuum model. Such a model might apply to a one- 
dimensional molecular crystal with a two-fold electronic degeneracy or near degeneracy. 
The existence of a double degeneracy is insured with sufficiently high spatial symmetry such 
as rotational symmetry axis three-fold or higher. The simplicity of the model considered here 
results in an analytic treatment of the wavef unctions, binding energies, and phase diagram. 

It was proposed lon g ag o that Jahn- Teller coupling plays an important role in both high- 



Tc superconductivity 



10| and fullerene superconductivity 



IjJ. More recently Jahn- Teller 



. Thus polarons of the kind described 



polarons have found application in the manganites 
here has relevance to several systems of current interest. 

On the basis of our model, we predict a new class of polarons in degenerate (or nearly 
degenerate) electron systems that couple to a bimodal distortion. We find that the model 
admits three types of stable vector soliton solutions. These soliton states consist of a localized 
electron with concomitant distortion of the elastic continuum - a polaron - having a spatial 
structure that involves a mixing of the degenerate (or nearly degenerate) orbital states. The 
electronic wavefunction of the solution is the direct analog of the envelope of the E-field 
in the vector solitons studied in non-linear birefringent optical fibers or the condensate 

n 

wavefunction in spinor Bose-Einstein condensation 

We observe that our study has applications to magnetic materials; one can reinterpret 
the two orbital components as those of a two-component spinor and apply our results to 
systems supporting magnetic polarons. 

The Model: We consider a single electron interacting with local, classical, bimodal dis- 
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tortions of the continuum. We take the Lagrangian to be the following 

C /^e -|- l~'ph ~l~ -^ep (l) 



with 



£e = - ^^9^^*) - \d^^^d.,^ - A^Va^ (2a) 



>Cp. = -Ei<^' (2b) 

i 

Cep = ^o</>o^^^ + ^1</>1*V3^ (2c) 

Ts is the usual Pauli matrix. 
We have used a two-component representation for our orbital basis with 



^ = (3) 




£e describes a doublet band of free carriers with splitting 2A. Cph describes two local 
classical distortions of the continuum with elastic constants kq and ki. We neglect the 
kinetic energy of the distortion, consistent with the adiabatic approximation. Cep gives the 
linear coupling of the free carriers to the local distortions: 0o preserves the local spatial 
symmetry, as a breathing mode would; 0i is a Jahn- Teller active mode that breaks the local 
symmetry and mixes the electronic states. We will refer to 0o as the Holstein mode, as it 
couples to the total electron density in the usual way of the Holstein model Atomic 
units are used where h = m = 1. 

The following equations of motion result from this Lagrangian: 

00 = ^(vl>tvl>) (4a) 



^(vl/tra^) (4b) 



We use Eq. |3]to eliminate distortions from the remaining equations of motion, yielding 
idtiJi - Ml + \dli'i + (^^iV^iT + r/|^2r)V^i = (5a) 
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idttp2 + + + {y\i^2? + v\i^i?)i^2 = 



(5b) 



with 



9o , 9i _ 2 , 2 
v = h — =7o + 7i 



(6a) 



9g 9i _ 2 2 



(6b) 



We observe that u is non-negative, while t] can change sign depending on the relative 
importance of the coupling to the Jahn- Teller mode 0i compared to the Holstein mode (pQ. 
In the context of the analogous system of equations in non-linear optics, u is referred to as 
the self-modulation coefficient, while rj is the cross-modulation coefficient. 

Consistent with Galilean invariance, we consider time-dependent solutions of the following 
form: 



E-A\ 



exp 



i{vx-{E+ -v^)t 



(7) 



A set of coupled nonlinear Schrodinger equations for the envelope functions rj results 



E-A 







(8a) 







(8b) 



Self-localized solutions exist only for E < and \E\ > A. In addition, known localized, 
stable solutions only exist for u and r] non- negative (l5l |. 
We define u"^ = {\E\ — A)/(|_E'| + A) and f3 = t]/u, giving 



r'l - ri + ri ( ri^ + /3r2 







(9a) 



(9b) 
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FIG. 1: Plot of type II polaron (a = 0.2 and uj = 0.6). Solid line is ri(^), and r2{£,) is dashed. 

with < /3 < 1 corresponding to the condition for known stable polaron solutions (see 
Fig. ED. 

Polaron Types: We describe three distinct types of stable polaron solutions to Eqs. IHi 
type I has equal amplitudes ri(^) = r2(^), corresponding to an equal orbital mixing that is 
constant over the spatial extent of the polaron. It results when uj = 1 and < /5 < 1. With 
a; = 1, the uncoupled electronic states must be exactly degenerate (A = 0). 



ri(0 = r2(0 = ^YT;^seche (10) 

The type II polaron has unequal amplitudes and results when oj ^ \. The value of 13 is 
determined by uj. The uncoupled electronic states are nearly degenerate. 

In the wave-daughter wave approximation, we take |ri| ^ \r2\- To leading order in 
ihe daughter wave amplitude, the following solution was determined by Yang and Benney 



ri(^) = V2 sech ^ (11a) 
r2(0 = a sech"^ ^ (lib) 

where a -C 1 and 

uj = ]^{^l^-l) (12) 
The type II polaron is unbound unless the coupling strengths satisfy 



7^ + 72 > v/24A 
5 



(13) 



Polaron type 


Binding Energy 


Type I 


24 '0 


Type II 




Type III 


24 '0 



TABLE I: Polaron binding energies. (Type II binding energy is given to 0{a).) 

Just as in the case of the impurity pseudo-Jahn- Teller effect, we observe that there is a 
critical coupling strength that must be exceeded for a distortion to develop [isf. 

In addition, Eq.^Jand the binding energy relation given in Table H] lead to the constraint 
that type II solutions exist only when 

P>-'-^ (14) 

This constraint consequently limits the range of u such that (\/2]~-|-T2v^ — 3)/6 < u < 1. 

The inequality of Eq. El and the binding energies of Table IJ lead to the boundaries of 
the region of stability for type II polarons in the model's parameter space. Equating the 
binding energies for type I and type II polarons gives the coexistence curve 



7? = ^J^| + 24A - 7o' (15) 

The type III polaron results in the case of (3 = 1 and u = 1. The equations in this 
case were studied first by Manakov 3] in another context. The condition of (3 = 1 implies 
that 7i must vanish. The special symmetry that results gives an infinite set of degenerate 
solutions, labeled by a continuous parameter < < 27r. This continuous internal degree of 
freedom allows one to qualitatively distinguish between the E ®hi polaron and the Holstein 
case. 

ri(^) = y2cos^ sech { (16a) 

r2(0 = V2sin^ sech ^ (16b) 

Fig. 121 gives the phase diagram for the three types of vector polarons. The upper region 
corresponding to 71 > 70 does not contain bound polarons. In the region 71 < 70, type I 
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FIG. 2: Plot of type III polaron (/3 = 1 and a; = 1) for 6* = f . Solid line is ri(^), and r2{0 is 
dashed. 

Yi 




FIG. 3: Plot of the polaron phase diagram. Type I polarons are stable in region I, and type II 
polarons are stable in region II. Type III polarons stably exist for 71 = (horizontal semi-axis). 
(We take the arbitrary case of A = 0.1.) 

polarons are stable throughout; however, type II polarons have a lower binding energy in 
the wedge-shaped region (light gray) . The type I polarons continue onto the locus of points 
7i = 0, becoming type III polarons. 

The one dimensional model studied has three types of stable two-component soliton 
solutions. These nonlinear excitations correspond to bimodal (pseudo) Jahn- Teller polarons, 
where degenerate (or nearly degenerate) orbital states are mixed under coupling to a localized 
distortion of the continuum. Such states might be produced in quasi-one dimensional systems 
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that are lightly doped by charge-transfer impurities. 

Type I and type III solutions have two important features in common: (1) the orbital 
mixing is constant in space, and (2) the wavefunction decay far from the polaron center has 
the same asymptotic form, going as exp(— We observe that type II solutions have a 
position-dependent mixing that results from the two components having different asymptotic 
forms, with the primary wavefunction going as exp(— 1^|) and the daughter wavefunction 
behaving as exp(— ti;|,^|). The amplitudes of the correlated distortions fall off twice as fast 
with distance. 

Finally we note that by interpreting the two states as the components of spin-| rather 
than those of an orbital doublet, the Lagrangian given in Eq. El also describes an electron 
interacting with local magnetization 0i. Thus these polarons can be thought of as new kinds 
of bound magnetic polarons in a system with an internal effective magnetic field of A which 
could be the result of exchange splitting or an externally applied field. Such states might 
have application to transport in doped magnetic semiconductors and other systems where 
there is interplay between electron-phonon coupling and magnetic exchange coupling. 
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